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Abstract. We have developed multi-dimensional constrained covariant density functional 



> 

!Tjr , theories (MDC-CDFT) for finite nuclei in which the shape degrees of freedom P\^ with even fi, 

f-^ ' e.g., /320, /522, /330, /?32, Pio, etc., can be described simultaneously. The functional can be one of 

VQ \ the following four forms: the meson exchange or point-coupling nucleon interactions combined 

with the non-linear or density-dependent couplings. For the pp channel, either the BCS approach 
or the Bogoliubov transformation is implemented. The MDC-CDFTs with the BCS approach for 
the pairing (in the following labelled as MDC-RMF models with RMF standing for "relativistic 
mean field") have been applied to investigate multi-dimensional potential energy surfaces and 
the non-axial octupole y32-correlations in TV = 150 isotones. In this contribution we present 
briefiy the formalism of MDC-RMF models and some results from these models. The potential 
energy surfaces with and without triaxial deformations are compared and it is found that the 
triaxiality plays an important role upon the second fission barriers of actinide nuclei. In the 
C^ , study of y32-correlations in Af = 150 isotones, it is found that, for ^""^Cf and ^^"Fm, 13^2 > 0.03 

and the energy is lowered by the /332 distortion by more than 300 keV; while for ■^'*®Cm and 
■^^^No, the pocket with respect to /332 is quite shallow. 



1. Introduction 

The shape of a nucleus can be described by the parametrization of the nuclear surface or the 
nucleon density distribution [1,2]. One of the mostly used parametrizations is the multipole 
expansion with coefficients /^a/^'s; see Fig. 1 [3] for a schematic show of some typical nuclear 
shapes. Many shape degrees of freedom are important not only for the ground states or small 
amplitude collective motions, but also for large amplitude collective motions such as fission. The 
readers are referred to Ref. [4] for an overview. Here we simply mention that both the non-axial 
and the refiection symmetries should be broken in the study of nuclear ground state properties 
and potential energy surfaces. 

We have developed multi-dimensional constrained covariant density functional theories 
(MDC-CDFT) by breaking the refiection and the axial symmetries simultaneously [3-7]. In 
these theories, the nuclear shape is assumed to be invariant under the reversion of x and y axes, 
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Figure 1. (Color online) A schematic show of some typical nuclear shapes. From left to right, 
the 1st row: (a) Sphere, (b) Prolate spheroid, (c) Oblate spheroid, (d) Hexadecapole shape, and 
the second row: (e) Triaxial ellipsoid, (f) Reflection symmetric octupole shape, (g) Tetrahedron, 
(h) Reflection asymmetric octupole shape with very large quadrupole deformation and large 
hexadecapole deformation. Taken from Ref. [3]. 



i.e., the intrinsic symmetry group is V4 and all shape degrees of freedom /3;^^ with even /U, e.g., 
/320) /322, /330) /332) Pao, ' ' 'j a™ included self-consistently. The covariant density functional can 
be one of the following four forms: the meson exchange or point-coupling nucleon interactions 
combined with the non-linear or density-dependent couplings. For the pp channel, either the BCS 
approach or the Bogoliubov transformation is implemented. The MDC-CDFT with the BCS 
approach for the pairing is named as MDC-RMF and that with the Bogoliubov transformation 
as MDC-RHB [4]. In this contribution, we will present the formalism for MDC-RMF models 
and some results of actinide nuclei. 

In Section 2, we give the formalism of MDC-RMF models briefly. The results for actinide 
nuclei are shown and discussed in Section 3. A summary is given in Section 4. 

2. Formalism of MDC-RMF models 

In this section, we briefly give the formalism of MDC-RMF models. More details can be found 
in Ref. [4]. The starting point of a RMF model with the non-linear point coupling interactions 
is the following Lagrangian [8-16], 



C = ij (i7/.5^ - Mb) V' - ^lin - >Cnl - >Cdcr - ^c 



(1) 
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(2) 



are the linear, non-linear, and derivative couplings and the Coulomb part, respectively. M^ is 
the nucleon mass, a^, ay, aTs-, CiTV-, f^s, Is, Iv-, ^s-, ^v-, ^TS-, and 5tv are coupling constants 
for different channels and e is the electric charge, ps, Pts-, Jv-, and jtv are the iso-scalar density, 
iso- vector density, iso-scalar current, and iso- vector current, respectively. 

Starting from the above Lagrangian and under several approximations, one can derive the 
Dirac equation for the nucleons. 



/iVi = {a • P + /3[Mb + S{v)] + V{v)} iPi = €iiPi, 



(3) 



where the potentials V{r) and S{r) are calculated from the densities. 

An axially deformed harmonic oscillator (ADHO) basis is adopted for solving the Dirac 
equation [3-7, 17-19]. Note that a RMF model with reflection asymmetry has been developed in 
a two-center HO basis [20] . The ADHO basis are defined as the eigen solutions of the Schrodinger 
equation with an ADHO potential. 
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where tOz and ujp are the oscillator frequencies along and perpendicular to the z axis, respectively. 
These basis are also eigen functions of the z component of the angular momentum jz with 
eigen values K = mi + nis- For any basis state ^^(ro"), the time reversal state is defined as 
<l'a(rc'") = T^a(ro"), where T = iUyK is the time reversal operator and K is the complex 
conjugation. Apparently we have Ka = —K^ and vr^ = vTq, with vr = ±1 being the parity. These 
basis form a complete set for expanding any two-component spinors. For a Dirac spinor with 
four components. 



V'i(ro-) 



(5) 



where the sum runs over all the possible combination of the quantum numbers a = 
{nz,nr,mi,ins} and /" and gf are the expansion coefficients. In practical calculations, one 
should truncate the basis in an effective way. 

In our model the nucleus is assumed to be symmetric under the V4 group, that is, for all the 
potentials and densities we can do the Fourier series expansion. 



f{p,V,z) = fo{p,z) 



1 1 

^ + H f^'^P' ^)^^ cos(2n99), 

2^ n=l ' ^ 



IT 



(6) 



where fo{p,z) and fnip,z) are real functions of p and z. The formalism for calculating the 
Fourier components of the potentials and densities can be found in Ref. [4]. 

Either the BCS approach or the Bogoliubov transformation has been implemented in our 
model to take into account the pairing effects. For the pairing force, we can use a delta force or 
a separable finite-range pairing force [21-23]. 

To obtain a potential energy surface, i.e., the energy of a nucleus as a function of deformation 
parameters, we make multi-dimensional constraint calculations which are equivalent to adding 
external potentials during the iteration. A modified linear constraint method was included in 
our MDC-RMF calculations [3-7]. 
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Figure 2. (Color online) Potential energy curve of ■^'''^Pu with various self-consistent symmetries 
imposed. The solid black curve represents the calculated fission path with V4 symmetry imposed, 
the red dashed curve that with the axial symmetry (AS) imposed, the green dotted curve that 
with the reflection symmetry (RS) imposed, the violet dot-dashed line that with both symmetries 
(AS & RS) imposed. The empirical inner (outer) barrier height is taken from Ref. [24] and 
denoted by the grey square (circle). The energy is normalized with respect to the binding 
energy of the ground state. The parameter set used is PC-PKl. Taken from Ref. [5]. 



The total energy of a nucleus is obtained by substituting the densities into the expectation 
value of the Hamiltonian. The center of mass correction Ec.m. can be calculated either 
phenomenologically or microscopically, depending on the effective interactions. The intrinsic 
multipole moments are calculated from the vector densities by 



Qx^ = J(frpv{ryYx^{n), 



(7) 



where Y\^{i}) is the spherical harmonics. 



3. Results and discussions 

In Ref. [5], one- (1-d), two- (2-d), and three-dimensional (3-d) constraint calculations were 
made for the actinide nucleus ^'^''Pu using the MDC-RMF models with the parameter set PC- 
PKl [25, 26]. In Fig. 2 we show the 1-d potential energy curve from an oblate shape with /32o 
about —0.2 to the fission configuration with /32o beyond 2.0 which are obtained from calculations 
with different self-consistent symmetries imposed: The axial (AS) or triaxial (TS) symmetries 
combined with reflection symmetric (RS) or asymmetric cases. The importance of the triaxial 
deformation on the inner barrier [27, 28] and that of the octupole deformation on the outer 
barrier are clearly seen: The triaxial deformation reduces the inner barrier height by more than 
2 MeV and results in a better agreement with the empirical value [24]; the RA shape is favored 
beyond the fission isomer and lowers very much the outer fission barrier. Besides these features, 
it was found for the first time that the outer barrier is also considerably lowered by about 1 MeV 
when the triaxial deformation is allowed. In addition, a better reproduction of the empirical 
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Figure 3. (Color online) Potential energy curves of ^^^U in the second barrier region calculated 
from MDC-RMF models with the NL3* [29,30], NL-Z2 [31], PC-PKl [25,26], DD-ME2 [32], 
and DD-PCl [33] parameter set, respectively. The reflection asymmetric shapes are allowed. 
The axial and non-axial symmetric results are denoted by black and red curves, respectively. 
Taken from Ref. [3]. 



barrier height can be seen for the outer barrier. It has been stressed that this feature can only 
be found when the axial and reflection symmetries are simultaneously broken [5]. 

In Ref. [5] , it was also examined the parameter dependence of the influence of triaxiality on 
the second flssion barrier and the lowering effect of the triaxiality on the second fission barrier 
was also observed when parameter sets other than PC-PKl are used. In Fig. 3 we show the 
comparison of the potential energy curves of ^^^U in the second barrier region calculated from 
the MDC-RMF models with different parameter sets and Lagrangian forms, including meson- 
exchange ones NL3* [29, 30], NL-Z2 [31], DD-ME2 [32], and point-coupling ones PC-PKl [25, 26] 
and DD-PCl [33]. The results with and without triaxiality are both presented. From the figure 
we can see that, although the absolute values of the barriers differ a lot among these results, for 
all five parameter sets the second barriers are lowered by the triaxiality. The largest effect is 
around 1 MeV for NL-Z2 and the smallest one is around 300 keV for DD-PCl. Thus we conclude 
that the lowering effects of the triaxiality on the second barriers are parameter independent in 
the RMF models. 

The self-consistent three-dimensional constraint calculations are very time-consuming. From 
the benchmark calculations for ^^''Pu, we learned many experiences about the important roles 
played by various shape degrees of freedom in different regions of the deformation space. For 
example, around the first fission barrier an actinide nucleus assumes triaxial and reflection 
symmetric shapes but around the second fission barrier both triaxial and octupole deformations 
are important [4]. These experiences are used in a systematic study of even-even actinide 
nuclei and the results were presented in Ref. [4]. Here in Fig. 4 we show the 1-d potential 
energy curves for even-even actinide nuclei from the fission isomer with /32o = 1-0 to the fission 
configuration with /32o = 1-7 calculated from the MDC-RMF model with the parameter set PC- 
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Figure 4. (Color online) Potential energy curves of even-even actinides nuclei in the second 
barrier regions calculated from the MDC-RMF model with the PC-PKl parameter set [25, 26]. 
The reflection asymmetric shapes are allowed. The axial and non-axial symmetric results are 
denoted by solid and dotted curves, respectively. The binding energies are normalized with 
respect to the ground states. The empirical values are taken from Ref. [24] and denoted by full 
circles. Taken from Ref. [3]. 



PKl [25, 26]. For comparison we present the results with and without the non-axial deformations. 
The reflection asymmetric shapes are allowed in both calculations. It was known that most of 
the outer barriers of the actinide nuclei are reflection asymmetric, which is consistent with 
the observed low-energy asymmetric fission fragments in this region. However, to estimate the 
fission barrier heights with an even higher accuracy, it is desirable to include also the effects of 
the non-axial deformations. From these potential energy curves it is clear that the non-axial 
deformations lower the outer barriers by around 0.5 to 1 MeV compared with the axial symmetric 
results. The empirical values of the fission barriers are also presented in the figure as full circles. 
Apparently, by including the triaxiality, the agreement between the theory and the experiment 
is improved. 

In Ref. [6] the non-axial reflection-asymmetric /332 shape in some transfermium nuclei with 
N = 150, namely ^'^'^Cm, ^'^^Cf, ^sopm^ ^nd ^^^No were investigated using the MDC-RMF 
model. In Fig. 5, we show potential energy curves for these N = 150 isotones. For ^^^Cm, 
the ground state deformation /332 = 0.020. The potential energy curve is rather flat around 
the minimum. We denote the energy difference between the ground state and the point with 
1^32 = by -Edepth which measures the energy gain with respect to the /332 distortion. For ^^^Cm, 
^depth is only 34 keV. For ^^Scf, ssopj^^ ^nd ^^^No, the minima locate at /332 = 0.037, 0.034, and 
0.025, respectively. The corresponding energy gain i^depth = 0.351, 0.328, and 0.104 MeV. As 
is discussed in Ref. [6] , the occurrence of the non-axial octupole ^32 correlations is mainly from 
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Figure 5. (Color online) The binding energy E (relative to the ground state) for N = 150 
isotones ^^^Cm (dashed line), ^^^Cf (solid line), ^^°Fm (long-dashed line), and ^^^No (dotted 
line) as a function of the non-axial octupole deformation parameter (3^2 ■ Taken from Ref. [6] . 



a pair of neutron orbitals [734]9/2 {I'ii^n) ^-nd [622]5/2 (1^(79/2) which are close to the neutron 
Fermi surface and a pair of proton orbitals [521]3/2 (vr/7/2) and [633]7/2 ('/rii3/2) which are close 
to the proton Fermi surface. 

It may be hard to conclude that these nuclei have static non-axial octupole deformations 
from these results because the potential energy curve is flat around the minimum and £^depth is 
small [6]. However, the present calculations at least indicate a strong y32-correlation in these 
nuclei. Both the non-axial octupole parameter (3^2 and the energy gain E'depth reach maximal 
values at ^^^Cf in the four nuclei along the A^ = 150 isotonic chain. This is consistent with the 
analysis given in Refs. [34, 35] and the experimental observation that in ^^*Cf, the 2~ state is 
the lowest among these nuclei [36] . 

4. Summary 

In this contribution we present the formalism and some applications of the multi-dimensional 
constrained relativistic mean field (MDC-RMF) models in which all shape degrees of freedom 
/3a^ deformations with even /x are allowed. The potential energy surfaces (curves) of ^^"^Pu and 
actinide nuclei with various symmetries are investigated. It is found that besides the octupole 
deformation, the triaxiality also plays an important role upon the second fission barriers. For 
most of even-even actinide nuclei, the triaxiality lowers the outer barrier by 0.5 ~ 1 MeV, 
accounting for about 10 ~ 20% of the barrier height. The non-axial reflection-asymmetric /332 
shape in some transfermium nuclei with A^ = 150, namely ^^^Cm, ^^^Cf, ^^"^Fm, and ^^^No are 
studied and rather strong non-axial octupole I32 effects have been found in ^^^Cf and ^^"^Fm 
which are both well deformed with large axial-quadrupole deformations, /32o ~ 0.3. We note 
that it is crucial to include the reflection asymmetric and non-axial shapes simultaneously for 
the study of potential energy surfaces and fission barriers of actinide nuclei and of nuclei in 
unknown mass regions such as superheavy nuclei. 



Acknowledgement 

This work has been supported by Major State Basic Research Development Program of China 
(Grant No. 2013CB834400), National Natural Science Foundation of China (Grant Nos. 
11121403, 11175252, 11120101005, 11211120152, and 11275248), Knowledge Innovation Project 
of Chinese Academy of Sciences (Grant Nos. KJCX2- EW-NOl and KJCX2-YW-N32). The 
results described in this paper are obtained on the ScGrid of Supercomputing Center, Computer 
Network Information Center of Chinese Academy of Sciences. 

References 



[1] 
[2] 

[3 



[4] 
[5] 
[6] 

[7] 



[8] 
[9] 

[lo; 

[11 

[12 

[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20^ 
[21 
[22 
[23 
[24 



[25 
[26' 
[27' 
[28; 

[29 

[3o; 

[31 
[32 
[33 
[34" 
[35' 
[36' 



Bohr A and Mottclson B R 1998 Nuclear Structure Vol I (World Scientific) 

Ring P and Schuck P 1980 The Nuclear Many-Body Problem ( Springer- Verlag Berlin, Heidelberg, and New 

York) 
Lu B N 2012 Multi- dimensional constrained relativistic mean field theory and the potential energy surfaces 

and fission barriers of actinide nuclei Ph.D. thesis. Institute of Theoretical Physics, Chinese Academy of 

Sciences (in Chinese) 
Lu B N, Zhao J, Zhao E G and Zhou S G 2013 (Preprint 1304.2513) 

Lu B N, Zhao E G and Zhou S G 2012 Phys. Rev. C 85 011301(R)-5 [Preprint 1110.6769) 
Zhao J, Lu B N, Zhao E G and Zhou S G 2012 Phys. Rev. C 86 057304-4 [Preprint 1209.6567) 
Lu B N, Zhao J, Zhao E G and Zhou S G 2012 EPJ Web of Conferences 38 05003, Proceedings of the 

Int. Conf. on Nuclear Structure and Related Topics (NSRT12), Dubna, Russia, 2-7 July 2012 [Preprint 

1303.0621) 
Serot B D and Walecka J D 1986 Adv. Nucl. Phys. 16 1-327 
Reinhard P G 1989 Rep. Prog. Phys. 52 439-514 
Ring P 1996 Prog. Part. Nucl. Phys. 37 193-263 

Vretenar D, Afanasjev A, Lalazissis G and Ring P 2005 Phys. Rep. 409 101-259 
Meng J, Toki H, Zhou S G, Zhang S Q, Long W H and Geng L S 2006 Prog. Part. Nucl. Phys. 57 470-563 

[Preprint nucl-th/0508020) 
Paar N, Vretenar D and Colo G 2007 Rep. Prog. Phys. 70 691-793 
Niksic T, Vretenar D and Ring P 2011 Prog. Part. Nucl. Phys. 66 519-548 
Prassa V, Niksic T, Lalazissis G A and Vretenar D 2012 Phys. Rev. C 86 024317-12 
Meng J, Peng J, Zhang S Q and Zhao P W 2013 Frontiers of Physics 8 55-79 
Lu B N, Zhao E G and Zhou S G 2011 Phys. Rev. C 84 014328-10 [Preprint 1104.4638) 
Gambhir Y K, Ring P and Thimet A 1990 Ann. Phys. 198 132-179 
Ring P, Gambhir Y K and Lalazissis G A 1997 Comput. Phys. Commun. 105 77-97 
Geng L S, Meng J and Toki H 2007 Chm. Phys. Lett. 24 1865-1868 
Tian Y, Ma Z and Ring P 2009 Phys. Lett. B 676 44-50 
Tian Y, Ma Z and Ring P 2009 Phys. Rev. C 80 024313-7 

Niksic T, Ring P, Vretenar D, Tian Y and Ma Z 2010 Phys. Rev. C 81 054318-12 
Capote R, Herman M, Oblozinsky P, Young P, Goriely S, Belgya T, Ignatyuk A, Koning A, Hilaire S, Plujko 

V, Avrigeanu M, Bersillon O, Chadwick M, Fukahori T, Ge Z, Han Y, Kailas S, Kopecky J, Maslov V, 

Reffo G, Sin M, Soukhovitskii E and Talou P 2009 Nucl. Data Sheets 110 3107-3214 
Zhao P W, Li Z P, Yao J M and Meng J 2010 Phys. Rev. C 82 054319-14 
Zhao P W, Song L S, Sun B, Geissel H and Meng J 2012 Phys. Rev. C 86 064324-6 
Abusara H, Afanasjev A V and Ring P 2010 Phys. Rev. C 82 044303-11 
Afanasjev A V 2013 Proceedings of Fifth International Conference on Fission and Properties of Neutron-rich 

Nuclei, November 4-10, 2012, Sanibel Island, USA [Preprint 1303.1206) 
Lalazissis G A, Konig J and Ring P 1997 Phys. Rev. C 55 540-543 
Lalazissis G A, Karatzikos S, Fossion R, Arteaga D P, Afanasjev A V and Ring P 2009 Phys. Lett. B 671 

36-41 
Bender M, Rutz K, Reinhard P G, Maruhn J A and Greiner W 1999 Phys. Rev. C 60 034304-20 
Lalazissis G A, Niksic T, Vretenar D and Ring P 2005 Phys. Rev. C 71 024312-10 
Niksic T, Vretenar D and Ring P 2008 Phys. Rev. C 78 034318-19 
Chen Y S, Sun Y and Gao Z C 2008 Phys. Rev. C 77 061305(R)-5 

Jolos R V, Malov L A, Shirikova N Y and Sushkov A V 2011 J. Phys. G: Nucl. Part. Phys. 38 115103 
Robinson A P, Khoo T L, Ahmad I, Tandel S K, Kondev F G, Nakatsukasa T, Seweryniak D, Asai M, Back 

B B, Carpenter M P, Chowdhury P, Davids C N, Eeckhaudt S, Greene J P, Greenlees P T, Gros S, Heinz 

A, Herzberg R D, Janssens R V F, Jones G D, Lauritsen T, Lister C J, Peterson D, Qian J, Tandel U S, 

Wang X and Zhu S 2008 Phys. Rev. C 78 034308-6 



